Transverse Presnel-Fizeau drag effects in strongly dispersive media. 
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A light beam normally incident upon an uniformly moving dielectric medium is in general sub- 
ject to bendings due to a transverse Fresnel-Fizeau light drag effect. In conventional dielectrics, 
the magnitude of this bending effect is very small and hard to detect. Yet, it can be dramatically 
enhanced in strongly dispersive media where slow group velocities in the m/s range have been re- 
cently observed taking advantage of the electromagnetically induced transparency (E1T) effect. In 
addition to the usual downstream drag that takes place for positive group velocities, we predict a 
significant anomalous upstream drag to occur for small and negative group velocities. Furthermore, 
for sufficiently fast speeds of the medium, higher order dispersion terms are found to play an im- 
portant role and to be responsible for peculiar effects such as light propagation along curved paths 
and the restoration of the spatial coherence of an incident noisy beam. The physics underlying this 
new class of slow-light effects is thoroughly discussed. 

PACS numbers: 42.50.Gy , 42.25.Bs 



I. INTRODUCTION 

A constant effort has always be devoted to the search 
for new effects and materials to control the propagation 
of light waves. Over the past few years, in particular, 
the use of quantum interference has led to an astonish- 
ing control of light waves propagating through specific 
classes of atomic and solid state media. These materials 
exhibit superior properties that cannot be found in con- 
ventional ones. Slow light propagation at group velocities 
as small as 1 m/s, e.g., has been observed in experiments 
with Bose- Einstein condensates of sodium atoms [l|, |2| , 
in hot rubidium vapors 0, H as well as in solid doped 
Pr:Y 2 Si05 crystals and Ruby @- Reversible stopping 
of a laser pulse Q in ultracold and hot alkali vapors HEjl 
as well as in Pr:Y2SiC>5 crystal has also been observed yg. 

Light stopping and slow-light propagation effects orig- 
inate from electromagnetically induced transparency 
(EIT). Such a widely discussed phenomenon arises from 
quantum interference and is characterized by a strong 
enhancement of the refractive index dispersion within 
a narrow frequency window around the medium reso- 
nance where absorption turns out to be largely quenched. 
The study of such a phenomenon goes back to the 
late seventies when non-absorbing resonances in atomic 
sodium have first been observed by Gozzini's group and 
later interpreted in terms of coherent population trap- 

pmg naojiiii. 

The interest for such a phenomenon has revived 0] 
over the past decade and has now become a rather top- 
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ical area of research [LJ, Il5| . much work being done 
on fundamental issues: high nonlinear coupling between 
weak fields and quantum entanglement of slow pho- 
tons 0, , entanglement of atomic ensembles |19| , 
quantum memories [2CJ and enhanced acousto-optical 
effects plj . just to mention a few. In particular, a 
strict analogy between slow-light in moving media and 
light propagating in curved space-times has been un- 
veiled and some of its consequences have been recently 
discussed HI 00 El 00] . 

Recently it has been also anticipated in [2^| that EIT 
media are ideal candidates for the observation of ex- 
tremely low and negative group velocities and (appar- 
ently) superluminal behaviour. In such media, the group 
velocity can in fact be readily tuned over a wide range 
of negative values directly by varying the coupling and 
probe detunings in a standard three-level A configura- 
tion. 

Although most slow light and negative group velocity 
experiments have dealt with the basic problem of a light 
pulse delay during its propagation across the dispersive 
medium, ultraslow positive or negative group velocities 
can have interesting consequences in many different sce- 
narios. In this paper we present a thorough investiga- 
tion of slow-light propagation through a moving medium. 
Specifically, we examine a configuration in which a highly 
dispersive EIT medium which moves with uniform veloc- 
ity and normally 1 to an incident light beam of finite spa- 



A configuration in which the medium is set to move parallel to 
the probe light-beam, leading to the more familiar longitudinal 
Fresnel-Fizeau effect, has been studied in l30l . Also in this con- 
figuration, the magnitude of the light-drag effect is predicted to 
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tial extent. The same geometry was adopted in the early 
seventies by Jones in his pioneering work on the trans- 
verse Fresnel-Fizeau light drag effect [^3, leading to 
the observation of a very small downstream bending, i.e., 
in the direction of motion, of a light ray. The use of a 
strongly dispersive medium supporting slow light, rather 
than the non-dispersive glassy material used by Jones, 
will not only allow for a remarkable enhancement of the 
drag effects but also for qualitatively new features. 

The possibility of having light propagating with small 
negative group velocities across the dragging medium is 
predicted to yield large upstream light bendings, i.e. in 
the direction opposite to that of the medium. The phe- 
nomenology of such an a anomalous Fresnel-Fizeau light 
drag effect, which has been the subject of an old contro- 
versy during the late seventies [33|, |34( , is here discussed 
in detail unwinding some of its controversial aspects. 

For sufficiently fast dragging speeds, the correct de- 
scription of the slow light transverse dragging effect re- 
quires that absorption dispersion and group velocity dis- 
persion be taken into account. As we shall see below, 
it turns out that these higher-order dispersion terms in- 
troduce new and peculiar features, such as propagation 
along curved light paths and the restoration of the spatial 
coherence of a noisy beam. All the numerical results pre- 
sented in this paper have been obtained using realistic 
parameters taken from slow-light experiments in ultra- 
cold atomic clouds [lj ■ Since the physics of the system 
is essentially determined by the electromagnetically in- 
duced transparency effect, the physical features here dis- 
cussed are however extremely general and hold through 
also for the recently prepared solid state EIT media. 

The paper is organized as follows. In sec. [H]we in- 
troduce the physical system and we present the model 
used for our predictions. The general theory of the trans- 
verse Fresnel-Fizeau light drag effect is presented in sec. 
IIIII These general results are specialized to the case of a 
strongly dispersive dressed medium driven into a A con- 
figuration by a resonant and non-resonant coupling beam 
respectively in sec. II VI and in sec. |V] in the former case, 
the group velocity is small and positive, so the transverse 
Fresnel drag occurs in the usual downstream direction; in 
the latter case, the negative group velocity is shown to 
give an anomalous upstream Fresnel drag. In sec I VII we 
proceed to discuss some effects which arise from the in- 
clusion of higher order dispersion terms. We carry out 
a complete analysis for two specific instances. In one 
case, the beam spectrum reshaping due to absorption 
and group velocity dispersion, which makes the average 
group velocity to bend during propagation, is seen to lead 
to curved light paths. In the other, the frequency disper- 
sion of absorption is seen to act as a filter for spatial 
fluctuations of the beam, whose space coherence proper- 
ties improve as it propagates through the medium. Fi- 



nally, a summary of the work is given in sec I VIII where 
conclusions are also drawn. 



II. THE MODEL AND GENERAL THEORY 

We consider a monochromatic probe light beam of fre- 
quency luq propagating along the z axis and normally in- 
cident upon a homogeneous dielectric medium uniformly 
moving along the x direction as shown in fig^ We here 
denote with L the thickness of the slab and with v its 
velocity v ^ c. The probe beam has a Gaussian profile 
centered in (xo,yo) so that at z — one has, 



E {x,y) = £ oe -[(— o) 2 +(^o) 2 ]/2^ ; 



(1) 



where (To is the beam waist. The corresponding Fourier 
transform Eo(k x ,k v ) of Eo(x,y) is then a Gaussian of 
width o~n . 




x,x 



Probe Beam 

FIG. 1: Scheme of the experimental set-up under considera- 
tion. 

In the following we shall always restrict our attention 
to the case of a weak probe beam, which allows us to 
apply linear response theory and describe the polariza- 
tion of the medium by means of a dielectric function 
e. Assuming for simplicity a non-magnetic (/x = 1) and 



be strongly enhanced in the slow light regime. 



isotropic e^.j = eSij medium, the scalar dielectric func- 
tion e completely characterizes the linear polarization of 
the medium in its rest frame Spatial locality of the 
dielectric polarization will also be assumed, i.e. e will be 
taken to depend on the frequency oj' but not on the wave- 
vector k'. Primed quantities will refer to the medium 
rest frame S', while the non- primed ones will refer to 
the laboratory frame E. The linearity of the optical re- 
sponse and the translational invariance of the system on 
the (x, y) plane allow us to decompose the incident field 
in its Fourier components in the (x 7 y) plane, to propa- 
gate each of them independently from the others, and to 
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finally reconstruct the transmitted beam profile by using 
an inverse Fourier transform. 

In the £' frame, the dispersion law inside the medium 
has the usual form: 

e{J)J 2 = c 2 {k' x 2 + k' y 2 + K 2 ). (2) 

Since the slab speed v <C c, the linearized form of the 
Lorentz transformations can be used, namely: 



k' 



u> — k x v 
k x - -xv 



(3) 
(4) 
(5) 



Notice that even though the medium does exhibits no 
spatial dispersion in its rest frame £', spatial dispersion 
arises however in the laboratory frame £ from the de- 
pendence of a/ onfcj 2 . By inserting the transformations 
(|3I5[) into the dispersion law J5J), one obtains the follow- 
ing expression: 



(uj - k x v) 2 



,(cj-k x v) = (k x -^-) 2 + k 2 y + k 2 z . (6) 



This equation represents the dispersion law in E and 
can be used to determine the propagation of the light 
beam across the slab. Knowing the frequency loq and the 
transverse components k x y of the wavevector allows us to 

obtain from @ the z component k z of the wave- vector 
inside the medium: 

kz ^ {k x j ky) — 



uj a - k x v 



Y e(uj - k x v) - (k x - - fc 2 , (7) 



c cr 
as well as in the external free space: 
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The profile of the transmitted beam can then be ob- 
tained by decomposing the field in its Fourier components 
at transverse wavevector (k x , k y ) and propagating each of 
them independently with the wave- vectors J7J and JSJ. 
Since the energy flux must be along the positive z axis, 



the roots with positive real part Re[fc 
be taken 3 . 



(in, out) -I 



> have to 



2 In the rest frame the dielectric polarization of a medium 
with frequency dispersion depends on the retarded values of the 
electric field at the same spatial position. As seen from the labo- 
ratory frame S, the polarization at a given point will thus depend 
on the electric field at different spatial positions, which means 
that a moving medium shows a spatial dispersion in S even if it 
is only temporally dispersive in the rest frame 

3 As the reflection amplitude is proportional to e — 1, the approxi- 
mation of neglecting interface reflections is a reasonable approx- 
imation in the case of an ultra cold atomic gas on which the 
present paper is focussed. On the other hand, a more complete 
theory including the possibility of multiple interface reflections 
is required in the case of a solid state medium. 



For each component, the amplitude at the position < 
z < L inside the slab is given by: 

E{k x ,k y -z) = e^ k *> k ^E (k Xl ky) (9) 

with a phase $: 

<S>{k x ,k y ;z) = kf n \k x ,k y )z (10) 

Notice that the wavevector k z m \k x ,k y ) as well as the 
phase <&{k Xl k y ] z) are generally complex quantities; their 
imaginary parts vanish only for non-absorbing, non- 
amplifying medium. Past the slab, the transmitted am- 
plitude at the position z > L is given by the same equa- 
tion © with the phase $ now given by: 

k y ; z) = ki m \k x ,k y ) L + k^ ut \k x ,k y ) (z-L). 

.(11) 

The spatial profile of the transmitted beam at any point z 
can then be obtained taking the inverse Fourier transform 
of ©, 

E( X , y, Z )= [ ^L^l. e i{k xX +k y y) e i*(k x , ky ;z) J, ). 

J 27T 

(12) 



III. THE TRANSVERSE FRESNEL-FIZEAU 
DRAG EFFECT 

Provided that the incident beam waist <7o is wide 
enough, only a very small window of wave- vectors (k x , k y ) 
around k x , y — is effectively relevant to the propagation 
dynamics and one can safely expand the phase (|llf> in 
powers of k x , y so that to the lowest order one has: 



$(k x ,k y ;z) = —y/e{uo) 



k x v 



e(w ) 



k x v de 
2 e(u>o) dto 



uj : 



L + ^(z-L). (13) 
c 



The position of the center (a: c , y c ) of the resulting wave 
packet at a given z can be obtained inserting the ex- 
pansion l|13|) into (|12fl and then invoke the so-called 
stationary-phase principle. This states that the integral 
in 112fl has its maximum value at those points (x c , y c ) 
for which constructive interference between the different 
Fourier components occurs, that is at those p_oints at 
which the phase of the integrand is stationary 



d 



dk x 



(Re[$(k X) k y ;z)] + k x 



kyVc) 



= 0. 



(14) 

Inserting into l|14f> the expression of the phase (|13fl and 
assuming the imaginary part q of e = e r + i to be negli- 
gible, one finds that in the laboratory frame S the beam 
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propagates inside the moving slab at a non- vanishing an- 
gle 6 with respect to the normal z direction given by: 



tan# 



w 



de r 



1 



dio 



_ ^ph 
c 



(15) 



Here v' gI and v' ph denote respectively the group and phase 
velocities in the medium rest frame E': 



A/e r (w ) + 



(16) 
(17) 



This deflection of the light beam can be interpreted 
as a transverse Fresnel-Fizeau drag effect, in which the 
beam of light is dragged by the transverse motion of the 
medium. After exiting from the rear surface of the slab, 
the beam again propagates along the normal direction. 
Its center, however, turns out to be laterally shifted by 
an amount, 



Ax = Lv 



1 



„2 



(18) 



v gi x/vgr z obtained from ijT§|) and (|2~l"1) agrees indeed with 

03- 

Another picture of this drag effect can be obtained 
by working in the rest frame E' j^]. Because of the 
aberration of light [3^ , the direction of the incident beam 
makes an angle 9' inc = —v/c with the normal 5 . In the 
rest frame E', the group and phase velocities inside the 
medium are parallel to each other and make an angle 
0' refT = —v/(cy/e) with the normal according to Snell's 
law. The modulus of the group velocity in E' is given by 
(|17fl . By Lorentz-transforming back the group velocity 
to the laboratory frame E, it is easy to check that the 
same deflection angle as in (|15|) is obtained. 

As one can see from the explicit expression l|18|) the 
magnitude of the transverse drag Ax is largest for a 
strongly dispersive medium in which e has a rapid fre- 
quency dependence and the group velocity v' gl . results 
much slower than the vacuum speed of light v' gI <C c. 

In this case, Galilean velocity composition laws can be 
safely applied and the result l|18l) simplifies to 



Ax = — , 

v 



(22) 



which yields a very intuitive interpretation for Ax as 
the displacement of the medium during the time inter- 
val At = L/v' gI taken by the light to travel across it. 



along a direction parallel to the medium velocity. This 
expression is in agreement with the ones derived by 
Player and by Rogers HEH^. 

While the group velocity direction inside the moving 
medium makes a finite angle with respect to the normal, 
the transverse x and y components of the phase velocity 
are always vanishing in the laboratory frame E because 
the light beam is normally incident on the slab and the 
transverse wave vector is conserved at the interface. The 
non-parallelism of the group and phase velocities arises 
then from the effective spatial dispersion acquired in the 
laboratory frame E by the moving medium 4 . Indeed, 
one can calculate the group velocity v gr = Vkw(k) at 
k x ,y = directly from the dispersion law JBJ to obtain 
for each component: 



v 1 



1 



1 dui 



= V 1- 



u gr,z 



(19) 

(20) 
(21) 



2^/1 du 



The non-vanishing value of v gr _ x is responsible for the 
transverse Fresnel-Fizeau drag. The value of tanf? = 



IV. THE TRANSVERSE FRESNEL-FIZEAU 
DRAG EFFECT IN A EIT MEDIUM 

The first experimental observation of the transverse 
Fresnel-Fizeau d rag effect was performed in the mid- 
1970's by Jones [31| using a rotating glass disk as mov- 
ing dielectric medium. In such a non-dispersive medium 
the phase and group velocities were of the order of c 
(c/v' gT w c/«p h « 1.5) and the limited rotation velocity 
of the disk at the laser spot (v ~ 2 x 10 4 cm/s) limited the 
lateral displacement to a distance of the order of 6 nm. 
However, a clever optical alignment technique allowed 
not only to observe the effect, but also to discriminate 
the validity of the result (|15fl from other possible expres- 
sions [3^. 

As remarked in the previous section, the use of a 
strongly dispersive medium allows for a significant en- 
hancement of the drag effect. Very slow group velocities 
can now be obtained in both atomic samples 0- El [1 
and solid-state media [a, Q by optically dressing a res- 
onant transition with a coherent coupling laser beam as 
shown in the A level scheme of fig|21 

Under the assumption of a weak probe beam, lin- 
ear response theory holds for the system dressed by the 



5 Notice that in the rest frame £', the beam waist appears as 
4 The assumed isotropy of the medium in its rest frame £' guar- uniformly translating along the negative x axis. This unfamiliar 

antees that group and phase velocities are always parallel in £'. feature however does not affect the argument. 
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FIG. 2: Scheme of the energy levels involved in the optical 
transitions. 



strong coupling beam: the resulting dielectric constant 
is to be interpreted as describing the linear response of 
the optically driven medium to the weak probe. Ob- 
viously, this dressed dielectric constant has a strongly 
nonlinear dependence on the coupling intensity so to 
account for all the nonlinear processes induced in the 
medium by the coupling field. In the simplest case of 
non-degenerate levels, the rest frame dielectric constant 
of our optically driven A configuration acquires the well- 
known form ei m m m , 
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|Oc| 2 



(23) 



where 0J e ( m ) and 7 e (m) are respectively the frequency and 
the linewidth of the excited e and metastable m states, 
where we have set uj g = 0. Here lo c is the frequency of 
the e «-> m coherent coupling beam 6 and fl c its Rabi 
frequency. The linewidth 7 m of the metastable state is 
much smaller than the linewidth "f e of the excited state. 
The / parameter quantifies the oscillator strength of the 
optical transition: for an ultra cold atomic gas at 
atomic densities n ~ 10 12 cm~ 3 , / is of the order of a few 
10~ 3 7 e . The background dielectric constant £oo takes 
into account the effect of all the other non-resonant tran- 
sitions and for an atomic gas can be taken to be 1 to a 
very good approximation. 

For a resonantly dressed medium, i.e. S c — uj c — (uj e — 
uj rn ) = 0, the dielectric function l|23|l in the neighborhood 
of the resonance \ui — w e \ <C "f e can be rewritten as: 



8nfi 

7e 



4Hg n 

7e ' 



(24) 



In the imaginary part of this expression shown in figEJi, 
it is easy to seize a narrow dip around to = u e in the 



6 Doppler shift of the coupling beam is avoided by choosing the 
direction of the coupling beam to be orthogonal to both the probe 
beam and the medium velocity. The waist of the coupling beam 
is taken as much larger than both the waist cro of the probe and 
the thickness L of the medium. 



i 

0,1 



(0 ' 


w 











(a>-&> e )/Y e 

FIG. 3: Resonantly dressed EIT medium at rest: plot of the 
imaginary (a) and real (b) parts of the dielectric function e 
and of the group velocity u gr (c). Medium parameters corre- 
spond to the case of an ultra cold 23 Na gas: j e ~ 2ty ■ 10 MHz, 
A e = 589 nm, / = 0.009 7,, ^ = 1, 7m = 10"%. The 
Rabi frequency of the resonant (8 C = 0) coupling beam is 
Q, c — 0.1 7 e . For this choice of parameters, the minimum 
(positive) group velocity is u gr w lm/s. 



otherwise wide absorption profile of the g — > e transi- 
tion. Provided r^/7e 3> 7 m , absorption at the center 
of the dip is strongly suppressed, yielding nearly perfect 
transparency for a resonant probe beam; this effect is the 
so-called electromagnetically induced transparency (EIT) 
effect lOj. The linewidth of the dip is T ~ ifl 2 c /j e and 
becomes strongly sub-natural (r <C 7e) for fl c <C 7e- The 
assumed inequality j m <C 7 e guarantees that a good level 
of transparency can be obtained simultaneously with a 
sub-natural linewidth of the dip. 

In the dip region, the real part (figEb) of the dielec- 
tric function l|23|l shows an extremely steep dispersion 
yielding the group velocity: 



' gr 
C 



in, 



2lTfuJ e 



(25) 



which can be orders of magnitude slower than the vacuum 
value c (figEt). This suggests that the importance of 
the transverse Fresnel-Fizeau drag effect l|18|l should be 
strongly enhanced in an EIT medium. 

To verify this prediction, a complete calculation of the 
profile of the transmitted probe beam can be numeri- 
cally performed by inserting the explicit expression of 
the dielectric function l|23|l into the dispersion law © 
and then performing the inverse Fourier transform i|12[l . 
For a probe exactly on resonance with the g — ► e transi- 
tion (u)q —w e ), the numerical result plotted in figEJi is in 
perfect agreement with the analytic prediction i|18f) when 
the value l|2*5)l for the group velocity is used. The near 
absence of absorption at the center of the dip guarantees 
that only a small fraction of the incident probe intensity 
is absorbed by the medium (figEJ)). 
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FIG. 4: Propagation of a light beam through a slowly moving 
(v = 0.01 m/s), resonantly dressed EIT medium. The optical 
parameters of the medium are the same as in fig|3] the ver- 
tical dashed lines correspond to the surfaces of the medium, 
whose thickness is taken as L = 93 /mi. The incident beam is 
resonant ujq = uj e and its waist is <to = 20 /mi. Downstream 
transverse Fresnel-Fizeau drag (a) and corresponding weak 
absorption of the beam (b). The x axis is oriented in the 
downstream direction. 



Group velocities as low as lm/s have been observed 
in ultra cold atomic gases. For the choice of parameters 
(w gr = lm/s and L « 100 /im) made in fig0J a medium 
velocity of the order of v = 0.01 m/s gives a lateral shift 
of the order of Ax « 1 /im, orders of magnitude larger 
than the one originally observed by Jones I22J • An 
even larger lateral shift Ax should be obtained by using 
a solid state material as dragging medium. Group veloc- 
ities as slow as v gr = 45 m/s have in fact been recently 
observed in a Pr doped Y 2 Si05 0] and Ruby @ crystals. 
Furthermore, the mechanical rigidity and the possibility 
of working at higher temperatures should allow one to 
study the drag effect on a thicker sample moving at a 
higher speed v. 



V. ANOMALOUS TRANSVERSE 
FRESNEL-FIZEAU DRAG EFFECT 

The discussion of the previous sections has focused on 
the most common case of media with normal dispersion. 
As Kramers-Kronig causality relations [3£j ensure that 



transverse Fresnel-Fizeau drag in the presence of anoma- 
lous dispersion, i.e. in the presence of negative group 
velocity. As negative group velocities in non-magnetic 
media are forbidden by Kramers-Kronig relations | 39| in 
all frequency regions where the medium is transparent 
and non-amplifying, it has been possible to observe neg- 
ative group velocities only in the presence of substantial 
absorption [40l l4ll li^] or in amplifying media ^| . No- 
tice that the negative group velocity effects which can 
be observed in left-handed media even in the absence 
of absorption arise from a completely different mecha- 
nism, that is from a simultaneously negative value of both 
the dielectric constant e and the magnetic susceptibility 
[i 0] . This class of effects are excluded from the present 
discussion, which is limited to non- magnetic (fx = 1) ma- 
terials. 

In all experiments performed up to now, negative 
group velocities have been demonstrated by observing 
that the pulse advances in time with respect to the same 
wave packet propagating in vacuum. This kind of appar- 
ently superluminal behavior obviously refers to the peak 
of the wave packet only, and not to the propagation ve- 
locity of information; as discussed in the review , the 
velocity of the front of a step-function signal can never 
exceed c. 
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FIG. 5: Plot of the imaginary part ei(uj) of the dielectric 
function (a), of its real part e r (w) (b), and of the group veloc- 
ity v gI (ui) (c) for a non-resonantly (S c — 10 7 e ) dressed EIT 
medium at rest. The other parameters are the same as in 

figEl 
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(26) 



for all frequency regions at which the medium is trans- 
parent and non-amplifying, the corresponding transverse 
Fresnel-Fizeau drag (|18l) results directed in the down- 
stream direction, as if the light were to be dragged by 
the moving medium. 

On the other hand, several papers during the 70's |33l 
|34| have discussed the possibility of having an upstream 



As we have recently anticipated in |38j, a significant 
upstream Fresnel-Fizeau drag effect should be observable 
in non-resonantly dressed EIT media, when the coupling 
beam is not exactly on resonance with the m — ► e transi- 
tion, but has a finite detuning |<5 C | = \u> m -\-u> c — w e \ > j e . 
This kind of media have in fact been shown to be good 
candidates for the observation of ultra slow and negative 
group velocities [2|j. In the case |<5 C | ^> j e , the dielectric 
function ilL'ol) in the neighborhood of ui — uu m + u c (figlSJ) 



7 



can be written in the Lorentzian form: 

47T/Q2 i 



e(u) = £oo j— 



4tt/ 
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w 2 - w - |(7 m + 



"?7e \ ■ 

^ (27) 

Largest absorptions occur at the Raman resonance with 
the two-photon transition from the ground g state to the 
metastable m state via the excited e state, i.e. at the 
two-photon resonance frequency 
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6, 



(28) 



The shift fl 2 / S c from the bare resonant frequency is due 
to the optical Stark effect induced by the coupling beam. 
The linewidth of the resonance line is the sum of the bare 
linewidth of the metastable m level plus a contribution 
which takes into account its decay via the excited e state: 



72 



n 2 



(29) 



For large coupling beam detunings \8 C \ ^S> £l c , the 
linewidth 72 results much smaller than the natural 
linewidth of the excited state, while the oscillator 
strength of the transition is itself weakened: 



_ n 2 



(30) 



The peak absorption (proportional to 72/72) does not 
vary for increasing values of |<5 c |/f2 c , at least as far as 
72 ^> 7m- On the other hand, the anomalous disper- 
sion at resonance u> — lu 2 is under the same conditions 
enhanced due to the narrower linewidth 72. The corre- 
sponding group velocity is given by 



1 gr 
C 



8TTfLU e n 2 



tit 



(31) 



which, in the limit 72 S> J m , is a factor r ) 2 /^5 2 slower in 
magnitude than the one obtained in the resonant 8 C = 
case considered in the previous section. This means 
that the use of a detuned coupling with \8 C \ > 7 e should 
enable one to observe significant upstream drags over an 
optical thickness L of the order of the absorption length, 
so that the transmitted probe intensity still remains an 
appreciable fraction of the incident one. 

In order to verify this expectation, we have again in- 
serted the explicit expression of the dielectric function 
(|23|l into the propagation equation JBJ and we have nu- 
merically performed the inverse Fourier transform so as 
to obtain the profile of the transmitted beam (figEt)- 
As expected, light bendings are found to occur in the 
upstream direction and the magnitude of the effect is 
in good agreement with the approximated analytical ex- 
pression l|18|) in which the imaginary part of e was ne- 
glected. We have also verified in figEb that the beam 



1 

0,8 
0,6 
0,4 
0.2 



(a) 


ill 
^^^^^^^^^^^^^^^^^^^^^^ 


- 




i , i , i 





-(b) ' 


i , | , | 






i , i , i 





z [um] 

FIG. 6: Propagation of a light beam through a slowly mov- 
ing (v = 0.01 m/s), non-resonantly dressed (<5 C = 10 j e ) EIT 
medium; probe frequency on resonance with the two-photon 
transition (ujq = The optical parameters of the medium 
are the same as in fig|S] the vertical dashed lines correspond 
to the surfaces of the medium, whose thickness is taken as 
L = 1.9 /jm. The incident beam waist is <7o = 20 /im. Anoma- 
lous upstream transverse Fresnel-Fizeau drag (a) and corre- 
sponding absorption of the beam (b). The x axis is oriented 
in the downstream direction. 



is not completely absorbed during propagation through 
the medium. Notice that for the same set of material 
parameters and the same experimental configuration as 
in the previous section, except for the coupling beam de- 
tuning <5 C , the group velocity |w gr | is now of the order 
of 0.01 m/s, i.e. a factor 100 slower than in the case of 
a resonant coupling. This explains the bigger deflection 
angle 9. 

Non-resonantly dressed EIT media are therefore good 
candidates for the experimental observation of significant 
upstream deflections and consequently large anomalous 
transverse Fresnel-Fizeau drags. Such an observation 
should finally resolve a controversy started in the sev- 
enties about the possibility of such an effect [33|, l34| . 
Furthermore, the anomalous upstream Fresnel-Fizeau 
drag could also provide an interesting way of demon- 
strating negative group velocities, in alternative to the 
usual |40l Ell l42l |43| measurement of the negative tem- 
poral delay of the pulse after its propagation across the 
medium. 



VI. HIGHER ORDER DISPERSION EFFECTS 

The effects discussed in the previous sections entirely 
rely on the strongly reduced value of the group velocity 
of light in EIT media. For slow enough medium veloci- 
ties, higher order dispersion effects such as the dispersion 
of absorption or the group velocity dispersion are indeed 
very small and can be hardly observed in the results plot- 
ted in figs01and|SJ On the other hand, these terms may 
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be no longer negligible for sufficiently large values of the 
medium velocity, regime in which the light propagation 
can exhibit qualitatively new features 7 . In the following 
of this section we shall discuss in detail two examples of 
such effects. 



A. Light propagation along a curved path 

In the present subsection, we shall discuss the effect of 
non-rectilinear light propagation due to simultaneously 
large dispersions of both absorption |4jj (proportional to 
dei/dui) and group velocity (proportional to d 2 e r /dui 2 ). 

For large enough values of the absorption disper- 
sion dei/duj, the absorption coefficient can have signifi- 
cant variations across the range of transverse k x vectors 
present in the incident beam. In this case, the center of 
mass of the k T wave- vector distribution 



kT(z) = 



Jdk x dky k x \E{k Xl ^y? 
J dk x dky \E(k x ,k y ; z)\ 2 



(32) 



moves from its initial value k x m (z = 0) = as the beam 
propagates through the medium. For a spatially wide 
Gaussian incident beam, we can limit the expansion of 
to the linear terms in k x<y and we find that the distribu- 
tion of transverse wave vectors keeps a gaussian shape, 
but its center of mass k x m (z) shifts to: 



dlm[k 



(in)-. 



dk x 



(33) 



For a negative value of dlm[ki ]/dk x , the k x > com- 
ponents result in fact less absorbed than the k x < ones 
so that the Gaussian spectrum shifts towards the k x m > 

region; vice versa for a positive value of <9Im[fci m ' ) ] / dk x . 

In the slow group velocity regime {v' <C c) , the depen- 
dence of the propagation wavevector yj on the transverse 
wave vector k XtV mainly comes from the Doppler effect 
combined with the strong frequency dispersion of the di- 
electric function, so that the transmission phase reads as: 



$(k x ,k y ;z) 



UJQZ 



y/e(ujQ - k x v). 



(34) 



By applying the same stationary-phase arguments used 
in sec IHII to the finite k x m (z) case, we find that after 



7 Even at the highest velocities here considered (v of the order 
of a few m/s), the use of the linearized form of the Lorentz 
transformations 13151 is well justified: as the magnitude of 
the transverse Doppler effect of the probe and coupling beam 
|Awj_| — i^o 1,2 / c2 i; 10 — 8 7e the contribution of the terms in v 2 
arising from the relativistic 7 factor is negligible as compared to 
the one coming from the dispersion of e. 



propagation for a distance z the spatial center of mass of 
the beam is located at 



Ax(z) 



(35) 



In physical terms, as the spectral center of mass k x m of 
the beam varies with z because of the filtering action of 
the absorption, the transverse Fresnel-Fizeau drag at a 
given position z has to be evaluated using the group ve- 
locity v' gr at the Doppler shifted frequency ujo — v k x m (z). 
In the presence of group velocity dispersion, the variation 
of k x m (z) with z implies that the light beam is no longer 
rectilinear, but acquires a finite curvature. Notice in par- 
ticular that the analytical expression for the angle 9(z) 
between the trajectory l|35|l and the normal z direction 



tan6>(z) 



dAx(z) 
dz 



dv' dk% 



v' gT (ui - vk x m ) [v' SI (u - vk™)] 



2 duj dz 



(36) 



now contains a term explicitly depending on the spectral 
center of mass shift dk x lrL /dz which was not present in 

03- 

A related effect was discussed in the time domain 
in |47j | where a small time-dependence of the group ve- 
locity for a light pulse propagating through a dispersive 
dielectric was predicted to show up as a consequence of 
the combined effect of the frequency dispersions of ab- 
sorption and of group velocity. 
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FIG. 7: Non-rectilinear light beam propagation through a 
moving (v = 4 m/s) EIT medium for a resonant coupling 
(S c — 0) and a slightly detuned probe (luq —cu e — 0.04 y e ); the 
incident beam waist is <to = 3 /im. (a) panel: curved beam 
path across the moving medium, (b) panel: intensity of the 
light beam at different depths in the medium. The optical 
parameters of the medium are the same as in fig0| the ver- 
tical dashed lines correspond to the surfaces of the medium, 
whose thickness is taken as L = 1.65 fim. 

Figs l7l8l report the result of numerical calculations for 
the specific case of a moving EIT medium when the cou- 
pling beam is resonant (S c — 0) and the probe beam is 
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k.v/y 



FIG. 8: Physical interpretation of the non-rectilinear light 
propagation of fig0 (a) panel: spatial Fourier transform of 
the incident (dashed) and transmitted (solid) beam profile; for 
each spectrum, the arrow indicates the position of the spec- 
tral center of mass. For the sake of clarity, the transmitted 
spectrum has been multiplied by 4. (b) and (c) panels: ab- 
sorption and inverse group velocity l/i> gr spectra as a function 
of the transverse wavevector k x . The corresponding Doppler- 
shifted frequency in the rest frame E' is uj' — ujq — k x v. Since 
ljq — oj e = 0.04 7 e , the resonance is found at k x v = 0.04 y e . 



tensity is of the order of 5% of the incident one (fig0)). 

It is also worth noticing that the curvature effect de- 
scribed in the present section follows from a reshaping 
of the beam in momentum space and hence is physically 
different from the ones discussed in [2^, Hj|, which in- 
stead originate from a non-uniform velocity field of the 
slow-light medium. 

B. Temporal and spatial coherence restoration 

If both probe and coupling beams are exactly on reso- 
nance (ojq — L0 e = 5 C = 0), both absorption dti/duj' and 
group velocity d 2 e r /du)' 2 dispersion vanish, so that the 
effects described in sec lVI Al do not take place. In the 
present subsection, we shall show how one can rather 
take advantage of the large value of cPei/dcu' 2 to im- 
prove the coherence level of a noisy incident probe beam. 
In sec I VI B Tl the case of temporal coherence restoration 
in a stationary EIT medium will be adressed, while in 
sec I VI B 2l we shall show how the same concepts can be 
applied in the case of a moving EIT medium to improve 
the level of spatial coherence. 

1. Temporal coherence restoration 



slightly detuned towards the blue (u>o — uj e = 0.04 7 e ). 
As one can verify by comparing the spectrum in figlS^ 
with the group velocity spectrum in figlSfc, most of the 
incident beam spectrum lies in the negative group ve- 
locity region (k x v < 0.02 j e ), so the beam is initially 
dragged in the upstream direction. As absorption is 
weaker for the k x > components for which the Doppler- 
shifted frequency lo' — ujq — k x v is closer to resonance 
(dha[kz ]/dk x < 0, see figEb), the negative k x < 
components are more rapidly quenched than the positive 
k x > components and the center of mass kx of the 
spectral distribution moves towards the positive k x > 
values. As one can see in figlHt, gf- -^r > in the region 

< k x v < 0.04 7 e of interest and therefore the curvature 
of the beam will be towards the downstream direction. 
After the first 1 /im of propagation, most of the k x spec- 
trum is found in the positive v' gr region (k x v > 0.02 j e ), 
so that the transverse Fresnel-Fizeau drag is from now 
on in the downstream direction (figGJi). 

As a consequence of the shift of the spectral center 
of mass k x m , the beams exits from the rear face of the 
medium at a small but finite angle with respect to the 
normal towards the downstream direction. Although this 
effect is hardly visible on the scale of figlEk., this small 
bending of the beam direction may have a significative 
effect on the subsequent rectilinear propagation of the 
beam in the free space. 

Unfortunately, this effect of non-rectilinear propaga- 
tion is associated to a rather severe absorption of the 
beam; for the specific case in figure, the transmitted in- 



Consider an incident probe beam of carrier frequency 
wo whose complex amplitude Eq (t) is assumed as fluctu- 
ating in time over a characteristic time scale r c : 

E{t) = E Q (t)e-^ ot . (37) 

Following a standard model [48(, the decay in time of 
the first-order coherence function g^(r) is taken to be 
Gaussian: 



9 {1) (r) = 



(E*(0)E (0)) 



exp(-r 2 /2r c 2 ). 



(38) 



The frequency spectrum |£?(uj)| 2 of the beam, being pro- 
portional to the Fourier transform of the coherence func- 
tion <?^(t), is also Gaussian with linewidth a c = 1/t c . 

In a one-dimensional geometry, the propagation of a 
pulse through a stationary dielectric medium is described 
by the usual Fresnel equation: 



(39) 



In the neighborhood of the resonance at uj e , the dielec- 
tric function of a slow light EIT medium can be approx- 
imately written as: 



e(w) ~ 1 + 



2c 

LU e V Rr 



(40) 



where w gr is the group velocity at resonance and the real 
quantity a is given by: 



d 2 6i{uj) 



dLU 2 



47r/7 e 



(41) 
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At the lowest order in u> — ui e , the real and imaginary 
parts of the wavevector k z (u>) can then be written as: 



0^ 

c 



1 + 



Wefgr 
^e) 2 



-{(jj - LU e ) 



(42) 
(43) 



In particular, notice how the imaginary part of k z is pro- 
portional to the square of w — uj e . Since after propagation 
over a distance L the amplitude of each frequency compo- 
nent is multiplied by a factor exp[ik z (u})L], the frequency 
spectrum after propagation keeps its Gaussian shape, but 
the frequency lincwidth is reduced to: 

a c {L) = °\ (44) 



1 



-L 



and the coherence time correspondingly increased to: 

T C (Z) = ■ 



1 



auj„ 



(45) 



Since the EIT medium does not provide amplification, a 
drawback of this filtering technique is that part of the in- 
cident intensity is lost during the line-narrowing process; 
the beam intensity after propagation through a distance 
L is in fact equal to: 



I(L) 



<Tc{L) 
(7c (0) 



7(0). 



(46) 



2. Spatial coherence restoration 

For a moving EIT medium and a strictly monochro- 
matic light beam at ujq, a similar effect occurs in the k x 
space. As discussed in full detail in sec[n] the propaga- 
tion of light is described in this case by equation Q . For 
resonant probe and coupling beams, the imaginary part 
of Imffcl" 1 '] corresponding to absorption is proportional 
to v 2 k1: 



lm[k z l 



4c 



1 2 l 2 

-v k, y 



(47) 



As the amplitude of the transverse spatial fluctuations is 
proportional to the amplitude of non-vanishing k x com- 
ponents, the spatial profile of the beam flattens as the 
beam propagates through the moving EIT medium. The 
faster the speed v of the medium, the more efficient the 
spatial coherence restoration process. 

Results of a numerical calculation for a spatially 
noisy incident beam propagating through a moving EIT 
medium are presented in figGD Notice how the fluctuation 
amplitude is strongly suppressed during propagation. At 
the end of the process, the losses in the total intensity 
amount to 60%. 

The overall shift of the beam which can be seen in 
the figure is due to the transverse Fresnel-Fizeau drag 
effect discussed in detail in sec II VI Since both the probe 
and the coupling are resonant, the shift is directed in the 
downstream direction. 
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FIG. 9: Transverse spatial coherence restoration during prop- 
agation across a moving (v — lm/s) EIT medium in a fully 
resonant regime (<5 C = luq — u e = 0, u>o — w e ). (a) panel: 
noisy incident beam profile. (b,c) panels: beam profile after 
propagation through the medium. Same medium parameters 
as in figl5H 



VII. CONCLUSIONS 

In the present paper we have given a comprehen- 
sive analysis of the transverse Fresnel-Fizeau drag ef- 
fect for light propagating across a uniform slab of mov- 
ing EIT medium. All calculations have been performed 
using realistic parameters taken from EIT experiments 
with ultra-cold atomic clouds. Since our results are es- 
sentially a consequence of electromagnetically induced 
transparency, they are extremely general and thus hold 
through also for the recently prepared solid state EIT 
media. Depending on the detuning of the probe and 
coupling beams with respect to the medium resonance, 
different regimes have been identified. 

In the presence of a slow and positive group velocity, 
the magnitude of the downstream Fresnel-Fizeau drag ef- 
fect is predicted to be significantly enhanced with respect 
to previous experiments. In the regime of negative group 
velocity, an anomalous upstream Fresnel-Fizeau drag has 
been predicted to occur. Not only would this help in 
solving a long-standing controversy on the observability 
of such an effect, but it could also provide an interesting 
alternative way for experimentally detecting a negative 
group velocity. 

For larger values of the velocity of the moving medium, 
higher order dispersion terms such as the dispersion of 
absorption and the dispersion of group velocity have been 
shown to play an important role in the phenomenology 
of light propagation. Depending on the specific choice 
of probe and coupling detuning, light propagation along 
curved paths can be observed, as well as the restoration 
of spatial coherence of a noisy beam. 

The extension of the present analysis to more compli- 
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catcd geometries including non-uniformly moving slow 
light media as well as probe pulses of finite duration will 
be the subject of future studies. 
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